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ABSTRACT 

Aims. We present results for the electron-impact excitation of highly-charged sulphur ions (S 8+ -S u+ ) obtained using the intermediate- 
coupling frame transformation ^-matrix approach. 

Methods. A detailed comparison of the target structure has been made for the four ions to assess the uncertainty on collision strengths 
from the target structure. Effective collision strengths (Ts) are presented at temperatures ranging from 2 x 10 2 (z + l) 2 K to 2 x 10 6 (z + 
l) 2 K (where z is the residual charge of ions). 

Results. Detailed comparisons for the Ys are made with the results of previous calculations for these ions, which will pose insight on 

the uncertainty in their usage by astrophysical and fusion modelling codes. 

Conclusions. 
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1. Introduction 

Many extreme ultraviolet (EUV) lines due to n = 2 — » 2 transi- 
tions of highly-charged sulphur ions have been recorded in so - 
lar observations (Thomas & Neupert [l994 Brown et al. 2008). 
Some transition lines show diagnostic potential for electro n den - 
sity, e.g. the emission lines of S XII (Keenan et al. 2002). 
Additionally, a few soft X-ray emission lines due to n — 3 — > 
2 transitions of sulphur ions have also been detected in solar ob- 
servations (Acton et al. 1985) and a Chandra Procyon observa- 
tion (Raassen et al. I2002T However, most of the excitation data 
adopted in astrophysical modelling for sulphur spectra are from 
the distorted-wave (DW) method, which is well known to under- 
estimate results for weaker transitions compared to those from 
an /^-matrix calculation. Only data from a few small 7?-matrix 
calculations is available for sulphur ions to-date. 

For S 8+ , Bhatia & Landi (2003a) reported an extensive exci- 
tation calculation (n=3) with the DW method, which is currently 
used by astrophysical modelling codes, e.g. chianti v6 (Dere et 
al. 120091) . /^-matrix data is available only for transitions within 
the ground configuration (Butler & Zeippen fl994l) . 

For S 9+ , DW excitation data for transi tions o f levels up to 
n — 3 were calculated by Bhatia & Landi (2003b), which were 
also incorporated into chianti v6. An i?- matrix calculation has 
been performed by Bell & Ramsbottom (2000), but only data 
for transitions among levels of n-2 complex and 2s 2 2p 2 3s con- 
figuration were reported. 

For S 10+ , the newest e xcitati on data are attributed to be the 
work of Landi & Bhatia d2003h . who have extended a previ- 
ous (n=3) DW calculation to include n=4 configurations (viz 



* These data are made available in the archives of APAP via 
http://www.apap-network.org OPEN-ADAS via http://open.adas.ac.uk 
as well as anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



2s 2 2p4Z', l'=s, p and d). Their results have been incorporated 
into astrophysical modelling codes, e.g. chianti v6. By using 
/^-matrix results available for some carbon-like ions (e.g. 2+ , 
Ne 4+ , Mg 6+ , Si 8+ and Ca 14+ ), Conlon et al. (fl992l) derived 
the electron excitation data for other carbon-like ions, includ- 
ing S I0+ , by interpolation. These resultant data are only valid 
for a range of electron temperatures approximately equal to 
logTmax + 0.8dex, where T max is the temperature of maximum 
fractional abundance in ionization equilibrium. However, inter- 
polation and/or extrapolation from 7?-matrix calculations was 
proved to be invalid at low temperatures because of the com- 
plexity of effective collision strengths along the iso-el ectronic 
sequence, as sh own i n Li-like (Liang & Badnel l 1201 1 1) . F-like 
(Witthoeft et al. l2007l). N e-like (Liang & Badnell |2010|f and Na- 
like (Liang et al. l2009bh iso-electronic sequences. An explicit 
/^-matrix calculation for this ion for transitions betwe en lev els 
of the n-2 complex was reported by Lennon & Burke dl994l) . 

For S 11+ , the / ^-mat rix calculations by Zhang et al. (1 19941) 
and Keenan et al. (2002) are still the main data source for mod- 
elling — the close-coupling expansion included the lowest 8 LS 
terms of the n-2 complex. However, there are potential prob- 
lems for B-like ions as demonstrated by Liang et al. (2009a) 
for Si 9+ , viz. effective coll ision s trengths of a previous /^-matrix 
calculation (Keenan et al. l2000l) do not converge to the correct 
high-temperature limit for a few strong transitions. 

Here, we report-on calculations for the electron-impact ex- 
citation of four iso-nuclear ions of sulphur (S 8+ — S 11+ ) which 
were made using the ICFT /^-matrix method. The remainder of 
this paper is organized as follows. In Sect. 2 and 3, we discuss 
details of the calculational method and pay particular attention 
to comparing our underlying atomic structure results with those 
of previous workers. The model for the scattering calculation 
is outlined in Sect. 4. The excitation results themselves are dis- 
cussed in Sect. 5 and we summarise in Sect. 6. 
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Table 2. Radial scaling parameters for S 8+ - S 11+ ions. 



Orbitals 


S s+ 




s iu + 


S u+ 


Is 


1.46525 


1.432 


1.39497 


1.38787 


2s 


1.25277 


1.339 


1.22509 


1.24802 


2p 


1.16771 


1.266 


1.15638 


1.17263 


3s 


1.21916 


1.232 


1.43868 


1.60433 


3p 


1.14234 


1.186 


1.26648 


1.36255 


3d 


1.16786 


1.281 


1.33814 


1.47452 


4s 


1.24628 


1.232 


1.36150 


1.47343 


4p 


1.15464 


1.189 


1.23441 


1.36462 


4d 


1.12290 


1.274 


1.30375 


1.53091 


4f 


1.13233 


1.289 


1.09803 


1.39692 



2. Structure: Level energy 

The target radial wavefunctions (Is - 4f) were obtained from 
autostructure (hereafter AS, Badnell 1986) using the Thomas- 
Femi-Dirac-Amaldi model potential. Relativistic effects were 
included perturbatively from the one-body Breit-Pauli opera- 
tor (viz. mass-velocity, spin-orbit and Darwin) without valence- 
electron two-body fine-structure operators. This is consistent 
with the operators included in the standard Breit-Pauli /^-matrix 
suite of codes. 

2.1. S IX 

Configuration interaction among 30 configurations (see Table[T} 
was included to describe the target used to calculate level en- 
ergies and weighted absorption oscillator strengths (gifij, for a 
given i <— j transition). The model potential radial scaling pa- 
rameters, A„i (n = 1 - 4; / e s, p, d, and f ), were obtained by a 
three-step optimization procedure. In the first step, the energy of 
the 2s A 2p v (x + y = 6) was minimized by varying the A\ s , /t_ s 
and A2 P scaling parameters. Then, the energies of the 2s 2 2p 3 3Z 
and 2s 2 2p 3 4/ configurations were minimized by varying the An 
and A41 scaling parameters, respectively. The resultant scaling 
parameters are listed in Table|2] 

The 92 lowest-lying fine-structure target levels of 2s A 2p- V (x+ 
y = 6), 2s 2 2p 3 3/ and 2s2p 4 3s ( 5 P and 3 P, only) configurations 
were used in the close-coupling expansion for the scattering cal- 
culation. The resultant level energies are compared with experi- 
mentally derived data from NIST v43 and previous calculations, 
see Table [3] The present AS level energies show an excellent 
agreement (less than 0. 5% except for 2s 2 2p 4 3 P) with those of 
Bhatia & Landi (2003a). This is due to the use of similar struc- 
ture codes and calculations. The main difference between the 
two is the n — 4 configurations that we inclu ded in our CI ex- 
pansion but were not used by Bhatia & Landi (2003a). 

When compared with N IST da ta and the MCHF collection^! 
(Tachiev & Froese Fischer, 2002), the present results agree to 
within 1% for all levels of the n=3 configurations. For levels of 
the n-2 configurations, the energy difference is about 2-5%. So, 
we performed a calculation with energy corrections to the diago- 
nal of the Hamiltonian matrix before diagonalization, for the 16 
lowest-lying levels, and iterated to convergence. For the miss- 
ing level (2s 2 2p 3 3s 5 S2) in the NIST compilation, we adopt the 
mean value of differences between our level energies and cor- 
responding NIST values of the same configuration. The resul- 
tant e-vectors and e-energies are used to calculate the oscillator 
strengths and archived energies. 

1 http://physics.nist.gov/PhysRefData/ASD/levels_form.html 

2 http://nlte.nist.gov/MCHF/ 



We notice that the MCHF data shows an excellent agreement 
with the data compiled in the NIST database. A non-relativistic 
multi-configuration Hartree-Fock (MCHF) approach was used 
by Tachiev & Froese Fischer (2002) to generate radial orbitals 
for subsequent use in diagonalizing a smaller scale Breit-Pauli 
Hamiltonian. A large set of configurations was used in the LS- 
coupling calculation, for example, configuration states were in- 
cluded up to n = 7. Orbitals sets were optimized separately for 
the initial and final states and no orthonormality was imposed 
between the two sets. Thus, we consider this MCHF data to be 
the theoretical reference work. 

2.2. S X 

Configuration interaction among 24 configurations (see Table[TJi 
was included to calculate the level energies and oscillator 
strengths between levels of the configurations 2s l 2p v (x +y — 5) 
and 2s 2 2p 2 3Z. Since the same configuration interaction and a 
simila r structu re code (superstructure) was used by Bhatia & 
Landi d2003bl) . we use their scaling parameters A„i in the present 
work for this ion. We note that there is severe interaction be- 
tween the 2s 2 2p 2 3Z and 2s2p 3 3Z' configurations. So some terms 
of the 2s2p 3 3s ( 6 S, 4 S and 4 D) and 2s2p 3 3p ( 6 P and 4 D) config- 
urations were included in the close-coupling expansion for the 
excitation calculation, as detailed in the next section. 

The calculated energies for the 84 lowest-lying levels are 
listed in Table|4]along with experimentally derived data from the 
NIST v4 compilation as well as other predictions. Although we 
have used the same con figuratio n expansion and same radial or- 
bitals as Bhatia & Landi ( 2003b), the energies do not quite match 
because we have omitted all two-body fine-structure operators 
so as to be consistent with our subsequent /^-matrix calculation. 
(If we include them for the 3 configurations of the ground com- 
plex then we reproduce their energies.) The difference is much 
smaller than the difference with the experimentally derived NIST 
energies. When compared with the NIST compilation^ and the 
MCHF collection El both agree to within 0.5% for all levels of 
the n — 3 configurations. For the levels of n — 2 complex, the 
difference is up to about 5%. So, we again iterated with energy 
corrections to the diagonal of Hamiltonian matrix before diago- 
nalization, for the 22 lowest-lying levels, to calculate oscillator 
strengths and archived en ergies. The data of MCHF collection 
(Tachiev & Froese Fischer |2002|) shows an excellent agreement 
again with the NIST data. 

2.3. S XI 

As shown in Table [TJ configuration interaction among 24 con- 
figurations has been taken into account to calculate the level en- 
ergies and oscillator strengths. The radial scaling parameters A„i 
were obtained by a three-step optimization procedure. In the first 
step, the energy of the 2s A 2p- v {x+y = 4) was minimized by vary- 
ing the /li s , ^2s an d ^2p scaling parameters. Then, the energies of 
the 2s 2 2p3/ and 2s 2 2p4Z configurations were minimized by vary- 
ing the An an d Mi scaling parameters, respectively. The resultant 
scaling parameters are listed in Table [2] 

The 254 lowest-lying fine-structure levels were used in the 
close-coupling expansion for the scattering calculation. They are 
compared with those data available from the NIST compilation 
and other predictions in Table [5] The present calculation shows 
a good agreement (1%) with those experimentally determined 
data in NIST database and the MCHF collection^! f or 2s2p 3 3 D, 
3 P and n = 3,4 levels. For other levels of the 2s2p 3 configura- 
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Table 1. Configurations included in the CC and CI expansions for S 8+ - S 11+ ions. 





CC configurations 


additional CI configuration 


C TV 


zs zp- 1 {x + y = o) 


Zs zp 4/, ZsZp 3s (except tor r ana r) 




zs zp 51 


ZsZp { jp, 3a, 4/J 




zszp Js ( r and r) 


Zp p,4)/ 






Zs Zp it 51' (x + y = Z) 


S X 


zs zp' {x + y = 5) 


ZsZp (3, 4)/ 




zs zp M 


Zs Zp 4/ 




zszp J s ( b, o and D) 


2p (3, 4}/ 




zszp jp { r ano r ) 




SXI 


2s*2p y (x + y = 4) 


2p J 3p('D, 'S) 




2s 2 2p3/, 2s2p 2 3/ 


2p 3 3p( 3 S, 13 P/D/F) 




2s 2 2p4/,2p 3 3s 


2s2p 2 4Z 




2p 3 3p (except for 'D, 'S) 


2p 3 4/ 




2p 3 3d (except for 3 S, 13 P/D/F) 




SXII 


2s x 2p y (x + y = 3) 


2p 2 3{4j/, 2s3s3/ 




2s 2 3(4)/, 2s2p3|4j/ 


2p3s3/, 2s3p 2 , 2s3d 2 



tion and those of the 2p 4 configuration, the present results are 
systematically higher than NIST data by 1-2%. The present AS 
result shows an excellen t agree ment (less than 0.5%) with the 
result of Landi & Bhatia d2003l) for all levels of the n-3 config- 
urations. However, both sets of results are systematically higher 
than the NIST data for the levels of n-2 complex, those of Landi 
& Bhatia (2003) more-so than the present which are within 2% 
(excluding the 5 S2). So, we perform an iterated energy correction 
calculation again for the 23 lowest-lying excited levels. 

The data of MCHF collection show better agreement again 
with the NIST data than other predictions. Unfortunately, there 
are no published papers to indicate the scale of calculations. 

2.4. S XII 

Configuration interaction among 32 configurations has been 
taken into account to calculate level energies and oscillator 
strengths, see Table [TJ The radial scaling parameters A„i were 
obtained by a three-step optimization procedure. In the first step, 
the energy of the 2s A 2p v (x + y = 3) configurations was mini- 
mized by varying the Au, z^s and fap scaling parameters. Then, 
the energies of the 2s 2 2p3/ and 2s 2 2p4/ configurations were min- 
imized by varying the A3; and A41 scaling parameters, respec- 
tively. The resultant scaling parameters are listed in Table [2] 

The 204 lowest-lying target levels were used in the close- 
coupling expansion for the scattering calculation. The present 
AS level energies are compared with those data available from 
the NIST compilation and other predictions, see Table|6] A good 
agreement (less than 1%) is obtained when compared with those 
experimentally determined in the NIST database for levels of the 
n=3 configurations. For levels of the n-2 complex, the differ- 
ence is slightly larger, but still within ~2%. Comparison with 
data in chianti v(0 demonstrates that the differences are within 
1% for almost all levels except for those of the n-2 complex. 
A good agreement is found when compared with calculations 
by Merkelis et al. d 19951) with many-body perturbation theory 
(MBP T) and all-order relativistic many-body theory by Nataraj 
et al. d2007l) . As done for other ions, energy corrections for the 
levels of 2s A 2p y (x + y-3) configurations have been included to 



3 Database comments denote that the gf compilation for S 11+ is from 
an unpublished calculation by Zhang et al. Additionally, the level ener- 
gies of the n=2 complex are the observed ones and not the theoretical 
values, as stated in this database. 



improve the accuracy of the oscillator strengths and archived en- 
ergies. 

3. Structure: Oscillator strengths 

A further test of our structure calculation is to compare weighted 
oscillator strengths gfy. In terms of the transition energy Ejj 
(Ryd) for the j —> i transition, the transition probability or 
Einstein's A-coefficient, Ap can be written as 

A i7 (au) = l^fElftj, (1) 

where a is the fine structure constant, and gt, gj are the statistical 
weight factors of the initial and final states, respectively. 

Figure [TJ shows such a comparison for the transitions into 
the five lowest-lying levels for the four iso-nuclear ions to assess 
the accuracy of the structure calculation. For S 8+ , about 8 6% of 
all available transitions in the work of Bhatia & Landi (2003a) 
show agreement to within 20%. When compared with the data 
from the MCHF collectioiJZl 61% of all available transitions 
agree to within 20%. For those transitions with larger differ- 
ences in the two cases, the data points are linked together by 
a solid line. For some transitions, the presen t results agree better 
with the results of Bhatia & Landi (2003a) than with the data 
from MCHF collection, while for others they agree better with 
the data from the MCHF method. Since correlation from much 
higher excited configuration has been taken into account in the 
data of MCHF collection, their gf- values are the best transition 
data so-far, as demonstrated by their level energies. The present 
results show a bett er agre ement with MCHF calculation (Tachie v 
& Froese Fischer, |2002|) than those of Bhatia & Landi (l2003al) . 
which indicates that we have a more accurate structure. 
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Table 3. Level energies (Ryd) of S from different calculations along with experimentally derived values from NIST v4. 



ID 


Specification 


NlSltH 


AS 


MCHF" 


BL03 6 


ID 


Specification 


NIST 


AS 


MCHF 


BL03 


1 


~ — — /I — TT^ 

2s 2 2p 4 J P 2 











47 


2s 2 2p 3 3d 5 D 4 




18.2653 


18.3135 


18.2392 


2 


2s 2 2p 4 3 P] 


0.07276 


0.0730 


0.0731 


0.0755 


48 


2s 2 2p 3 3p 'Pi 




18.2871 


18.3379 


18.2616 


3 


2s 2 2p 4 3 P 


0.097032 


0.0981 


0.0968 


0.1005 


49 


2s 2 2p 3 3p 3 P 2 




18.3578 


18.4164 


18.3383 


4 


2s 2 2p 4 'D 2 


0.531213 


0.5601 


0.5353 


0.5625 


50 


2s 2 2p 3 3p 3 Pi 




18.3728 


18.4217 


18.3583 


5 


2s 2 2p 4 'So 


1.1181 


1.1229 


1.1177 


1.1318 


51 


2s 2 2p 3 3p 3 P 




18.3808 


18.4251 


18.3685 


6 


2s2p 5 3 P 2 


4.05502 


4.1504 


4.0671 


4.1505 


52 


2s 2 2p 3 3d 3 D 2 


18.5464 


18.5248 


18.5519 


18.5046 


7 


2s2p 5 3 P! 


4.11888 


4.2172 


4.1312 


4.2165 


53 


2s 2 2p 3 3d 3 D! 


18.5475 


18.5265 


18.5542 


18.5065 


8 


2s2p 5 3 P 


4.15438 


4.2528 


4.1669 


4.2530 


54 


2s 2 2p 3 3d 3 D 3 


18.5522 


18.5321 


18.5585 


18.5119 


9 


2s2p 5 'P! 


5.61407 


5.8016 


5.6353 


5.7982 


55 


2s 2 2p 3 3p 'D 2 




18.6042 


18.5960 


18.5981 


10 


2p 6 'So 


9.470060 


9.7908 


9.5003 


9.7967 


56 


2s 2 2p 3 3p 'So 




18.9893 


18.9748 


18.9813 


11 


2s 2 2p 3 3s 5 S 2 




15.9343 


16.0032 


15.9188 


57 


2s 2 2p 3 3d 3 F 2 




19.0114 


19.0399 


18.9807 


12 


2s 2 2p 3 3s 3 Sj 


16.2493 


16.2047 


16.2520 


16.1900 


58 


2s 2 2p 3 3d 3 F 3 




19.0287 


19.0537 


18.9982 


13 


2s 2 2p 3 3s 3 D l 


16.8126 


16.7815 


16.8207 


16.7656 


59 


2s 2 2p 3 3d 'So 




19.0460 


19.0621 


19.0154 


14 


2s 2 2p 3 3s 3 D 2 


16.8144 


16.7845 


16.8225 


16.7625 


60 


2s 2 2p 3 3d 3 F 4 




19.0493 


19.0701 


19.0191 


15 


2s 2 2p 3 3s 3 D 3 


16.8208 


16.7950 


16.8289 


16.7761 


61 


2s 2 2p 3 3d 3 G 3 




19.1082 


19.1308 


19.0799 


16 


2s 2 2p 3 3s 'D 2 


16.9386 


16.9211 


16.9480 


16.9024 


62 


2s 2 2p 3 3d 3 G 4 




19.1145 


19.1344 


19.0864 


17 


2s 2 2p 3 3p 5 P! 




16.9639 


17.0182 


16.9377 


63 


2s 2 2p 3 3d 3 G 5 




19.1218 


19.1377 


19.0939 


18 


2s 2 2p 3 3p 5 P 2 




16.9696 


17.0242 


16.9438 


64 


2s 2 2p 3 3d 'G 4 




19.1466 


19.1597 


19.1210 


19 


2s 2 2p 3 3p 5 P 3 




16.9795 


17.0352 


16.9542 


65 


2s 2 2p 3 3d 3 D, 


19.1852 


19.1769 


19.2311 


19.1507 


20 


2s 2 2p 3 3s 3 P 
2s 2 2p 3 3s 3 Pi 




17.1515 


17.2265 


17.1317 


66 


2s 2 2p 3 3d 3 D 2 


19.2106 


19.2066 


19.2224 


19.1816 


21 




17.1577 


17.2315 


17.1378 


67 


2s 2 2p 3 3d 3 D 3 


19.2117 


19.2152 


19.2257 


19.1912 


22 


2s 2 2p 3 3s 3 P 2 


17.2396 


17.1736 


17.2447 


17.1536 


68 


2s 2 2p 3 3d 'P! 


19.2177 


19.2191 


19.1974 


19.1937 


23 


2s 2 2p 3 3p 3 Pi 




17.2128 


17.2521 


17.1940 


69 


2s 2 2p 3 3d 'D 2 


19.2954 


19.2740 


19.3080 


19.2480 


24 


2s 2 2p 3 3p 3 P 2 




17.2171 


17.2555 


17.1983 


70 


2s 2 2p 3 3d 3 P 2 


19.2865 


19.2876 


19.3004 


19.2630 


25 


2s 2 2p 3 3p 3 P 




17.2184 


17.2585 


17.1999 


71 


2s 2 2p 3 3d 3 P! 


19.3128 


19.3153 


19.3249 


19.2908 


26 


2s 2 2p 3 3s 'Pi 


17.3533 


17.2987 


17.3625 


17.2788 


72 


2s 2 2p 3 3d 3 P 




19.3178 


19.3283 


19.2929 


27 


2s 2 2p 3 3p 'Pi 




17.6866 


17.7191 


17.6583 


73 


2s 2 2p 3 3d 3 Si 


19.3693 


19.3521 


19.3829 


19.3291 


28 


2s 2 2p 3 3p 3 D 2 




17.7395 


17.7706 


17.7221 


74 


2s 2 2p 3 3d 'F 3 


19.4529 


19.4622 


19.4678 


19.4390 


29 


2s 2 2p 3 3p 3 D, 




17.7499 


17.7808 


17.7109 


75 


2s2p 4 3s 5 P 3 




19.4781 






30 


2s 2 2p 3 3p 3 D 3 




17.7588 


17.7866 


17.7305 


76 


2s 2 2p 3 3d 3 F 4 




19.4928 


19.5456 


19.4632 


31 


2s 2 2p 3 3p 3 F 2 




17.8168 


17.8469 


17.7910 


77 


2s 2 2p 3 3d 3 P 


19.5264 


19.4967 


19.5583 


19.4673 


32 


2s 2 2p 3 3p 3 F 3 




17.8289 


17.8575 


17.8032 


78 


2s 2 2p 3 3d 3 F 2 




19.5004 


19.5554 


19.4708 


33 


2s 2 2p 3 3p 3 F 4 




17.8422 


17.8683 


17.8169 


79 


2s 2 2p 3 3d 3 F 3 




19.5060 


19.5542 


19.4798 


34 


2s 2 2p 3 3p 'F 3 




17.8742 


17.9008 


17.8518 


80 


2s 2 2p 3 3d 3 Pi 


19.5449 


19.5118 


19.5705 


19.4828 


35 


2s 2 2p 3 3p 3 S! 




18.0829 


18.2133 


18.0569 


81 


2s 2 2p 3 3d 3 P 2 


19.5613 


19.5296 


19.5851 


19.5004 


36 


2s 2 2p 3 3p 3 P 




18.1322 


18.1492 


18.1138 


82 


2s2p 4 3s 5 P 2 




19.5299 






37 


2s 2 2p 3 3p 3 P 2 




18.1596 


18.1624 


18.1468 


83 


2s2p 4 3s 5 Pi 




19.5629 






38 


2s 2 2p 3 3p 3 Pi 




18.1770 


18.1262 


18.1570 


84 


2s 2 2p 3 3d 3 D 2 


19.6340 


19.6039 


19.6485 


19.5786 


39 


2s 2p 3p Dj 




18.1936 


18.2533 


18.1659 


85 


2s 2p 3d D 3 


19.6508 


19.6167 


19.6642 


19.5907 


40 


2s 2 2p 3 3p 3 D 2 




18.1958 


18.2505 


18.2390 


86 


2s 2 2p 3 3d 3 Di 


19.6493 


19.6183 


19.6651 


19.5936 


41 


2s 2 2p 3 3p 3 D 3 




18.2102 


18.2710 


18.1828 


87 


2s 2 2p 3 3d 'D 2 


19.7085 


19.6963 


19.7220 


19.6733 


42 


2s 2 2p 3 3p 'D 2 




18.2561 


18.2975 


18.1694 


88 


2s 2 2p 3 3d 'F 3 


19.7429 


19.7241 


19.7557 


19.7042 


43 


2s 2 2p 3 3d 5 D 




18.2617 


18.3130 


18.2352 


89 


2s2p 4 3s 3 P 2 




19.8880 






44 


2s 2 2p 3 3d 5 D! 




18.2620 


18.3131 


18.2356 


90 


2s2p 4 3s 3 Pi 




19.9453 






45 


2s 2 2p 3 3d 5 D 2 




18.2627 


18.3131 


18.2364 


91 


2s 2 2p 3 3d 'Pi 


19.9588 


19.9473 


19.9727 


19.9269 


46 


2s 2 2p 3 3d 5 D 3 




18.2638 


18.3131 


18.2375 


92 


2s2p 4 3s 3 P 




19.9718 







a Refers to data from the MCHF collection^!. 



* Refers to the work of Bhatia & Landi ( 12003 st) . 



For S , 82% of transitions agree to with in 20% for the 
present AS results and those of Bhatia & Landi (2003b). Recall, 
we omitted two-body fine-structure but have iterated to the ob- 
served energies, for leve ls of the ground complex, compared to 
Bhatia & Landi (2003b). We have also performed calculations 
with/without the two-body fine-structure and level energy cor- 
rections to study the effect of the two, see Table [7] It appears 
that the energy corrections play an more important on the resul- 
tant gf '-value for these weak transitions. When compared with 

the data from the MCHF collection^, about 64% of all avail- 
able transitions show agreement to within 20%. For those transi- 
tions with larger differences, the data points are linked together 



as done for S 8+ . We also notice that the level label for 2s 2 2p 2 3s 
2 D 5/2 (21 th in the Table© and 2 D 3/2 (22 th ) in MCHF collec- 
tion should be exchanged because a good agreement between the 
MCHF calculation and the other two predictions can be obtained 
after such a procedure for all transitions to the five lowest-lying 
levels from the two levels, e.g. the 22— > 3 and 2 1 — >3 transitions 
marked in Fig.Q] 

For S 10+ , most (67%) transitions are in agreement to 
within 20% f or the present AS results and those of Landi & 
Bhatia d2003l) . When compared with calculation from MCHF 
method, the percentage is about 90% of available transition data. 
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Table 4. Level energies (Ryd) of S from different calculations along with experimentally derived values from NIST v4. 



ID 


Specification 




AS 


BL03 fl 


MCHF* 


ID 


Specification 


NIST 


AS 


BL03 


MCHF 


1 


2s 2 2p 3 4 S 3/2 











43 


2s 2 2p 2 3p 2 P 3/2 




21.0732 


21.0731 


21.0498 


2 


2s 2 2p 3 2 D 3/2 


0.751270 


0.7929 


0.7963 


0.7568 


44 


2s 2 2p 2 3d 4 F 3/2 




21.1123 


21.1123 


21.1308 


3 


2s 2 2p 3 2 D 5/2 


0.761773 


0.8101 


0.8079 


0.7681 


45 


2s 2 2p 2 3d 4 F J/2 




21.1338 


21.1337 


21.1511 


4 


2s 2 2p 3 2 P 1/2 
2s 2 2p 3 2 P 3/2 


1.15708 


1.1932 


1.1944 


1.1597 


46 


2s 2 2p 2 3d 4 F 7/2 




21.1661 


21.1660 


21.1818 


5 


1.17375 


1.2131 


1.2119 


1.1765 


47 


2s 2 2p 2 3d 4 F 9/ 2 




21.2102 


21.2101 


21.2225 


6 


2s2p 4 4 P 5/2 


3.44876 


3.4966 


3.4954 


3.4608 


48 


2s2p 3 3s 4 S 3/2 




21.2346 




21.2768 


7 


2s2p 4 4 P 3/2 


3.51168 


3.5593 


3.5608 


3.5239 


49 


2s 2 2p 2 3d 2 P 3/2 


21.2431 


21.2356 


21.2355 


21.2470 


8 


2s2p 4 4 P 1/2 


3.54376 


3.5933 


3.5942 


3.5557 


50 


2s 2 2p 2 3d 4 D, /2 




21.2421 


21.2420 


21.2611 


9 


2s2p 4 2 D 3/2 


4.74518 


4.8730 


4.8739 


4.7647 


51 


2s 2 2p 2 3d 4 D 5/2 




21.2650 


21.2649 


21.2804 


10 


2s2p 4 2 D 5/2 


4.74646 


4.8771 


4.8763 


4.7660 


52 


2s 2 2p 2 3d 4 D 3/2 




21.2826 


21.2825 


21.2943 


11 


2s2p 4 2 S 1/2 


5.54765 


5.6861 


5.6864 


5.5645 


53 


2s 2 2p 2 3d 4 D 7/2 




21.2877 


21.2877 


21.3020 


12 


2s2p 4 2 P 3/2 


5.80384 


5.9722 


5.9715 


5.8300 


54 


2s 2 2p 2 3p 2 P 3/2 
2s 2 2p 2 3p 2 Pi /2 




21.3270 


21.3269 




13 


2s2p 4 2 P 1/2 


5.88369 


6.0544 


6.0545 


5.9101 


55 




21.3303 


21.3302 




14 


2p 5 2 P 3/2 


9.03293 


9.2852 


9.2860 


9.0635 


56 


2s 2 2p 2 3d 2 P 1/2 


21.3346 


21.3351 


21.3350 


21.3404 


15 


2p 5 2 Pi/ 2 


9.134286 


9.3930 


9.3914 


9.1650 


57 


2s 2 2p 2 3d 2 F 5/2 


21.3509 


21.3598 


21.3597 


21.3576 


16 


2s 2 2p 2 3s 4 P 1/2 


19.0223 


18.9961 


18.9959 


19.0238 


58 


2s 2 2p 2 3d 4 P 5/2 


21.4199 


21.4117 


21.4115 


21.4157 


17 


2s 2 2p 2 3s 4 P 3/2 


19.0674 


19.0380 


19.0379 


19.0660 


59 


2s 2 2p 2 3d 4 P 3/2 


21.4441 


21.4359 


21.4358 


21.4399 


18 


2s 2 2p 2 3s 4 P 5/ 2 


19.1224 


19.0981 


19.0979 


19.1226 


60 


2s 2 2p 2 3d 2 F 7/2 


21.4283 


21.4375 


21.4374 


21.4341 


19 


2s 2 2p 2 3s 2 P 1/2 


19.2560 


19.2446 


19.2445 


19.2535 


61 


2s 2 2p 2 3d 4 P 1/2 


21.4491 


21.4484 


21.4482 


21.4518 


20 


2s 2 2p 2 3s 2 P 3/2 
2s 2 2p 2 3s 2 D 5/2 
2s 2 2p 2 3s 2 D 3/2 
2s 2 2p 2 3p 2 S 1/2 
2s 2 2p 2 3p 4 D 1/2 
2s 2 2p 2 3p 4 D 3/2 
2s 2 2p 2 3p 4 D 5/2 
2s 2 2p 2 3p 4 Pi/2 
2s 2 2p 2 3p 4 D 7/2 
2s 2 2p 2 3p 4 P 3/2 
2s 2 2p 2 3p 4 P 5/2 


19.3234 


19.3163 


19.3161 


19.3221 


62 


2s2p 3 3p 6 P 3/2 




21.5973 




21.7238 


21 


19.6846 


19.6882 


19.6881 


19.6950 


63 


2s2p 3 3p 6 P 5/2 




21.6051 




21.7299 


22 


19.6768 


19.6907 


19.6905 


19.6924 


64 


2s2p 3 3p 6 P 7/2 




21.6169 




21.7420 


23 




19.8772 


19.8771 


19.8971 


65 


2s 2 2p 2 3d 2 D 3/2 


21.6763 


21.6959 


21.6957 


21.6829 


24 




19.9500 


19.9499 


19.9770 


66 


2s 2 2p 2 3d 2 D 5/2 


21.6872 


21.7147 


21.7145 


21.6984 


25 




19.9719 


19.9718 


19.9993 


67 


2s 2 2p 2 3d 2 G 7/2 




21.8602 


21.8601 


21.8540 


26 




20.0119 


20.0118 


20.0389 


68 


2s 2 2p 2 3d 2 G 9/2 




21.8693 


21.8692 


21.8617 


27 




20.0652 


20.0651 


20.0927 


69 


2s 2 2p 2 3d 2 D 3/2 


21.9401 


21.9623 


21.9622 


21.9494 


28 




20.0659 


20.0658 


20.0902 


70 


2s 2 2p 2 3d 2 D 5/2 


21.9533 


21.9702 


21.9701 


21.9512 


29 




20.0742 


20.0741 


20.1002 


71 


2s2p 3 3p 4 P 5/2 


22.0342 


21.9714 




22.0584 


30 




20.1104 


20.1103 


20.1344 


72 


2s2p 3 3p 4 P 3/2 




21.9756 




22.0633 


31 


2s 2 2p 2 3p 2 D 3/2 




20.2002 


20.2001 


20.2102 


73 


2s2p 3 3p 4 P 1/2 




21.9795 




22.0670 


32 


2s 2 2p 2 3s 2 S 1/2 




20.2363 


20.2361 


20.2976 


74 


2s 2 2p 2 3d 2 F 7/2 


21.9531 


21.9922 


21.9921 


21.9623 


33 


2s 2 2p 2 3p 4 S 3/2 




20.2660 


20.2658 


20.2949 


75 


2s 2 2p 2 3d 2 F J/2 


21.9848 


22.0203 


22.0201 


21.9959 


34 


2s 2 2p 2 3p 2 D 5/2 
2s 2 2p 2 3p 2 P I/2 
2s 2 2p 2 3p 2 P 3/2 




20.2749 


20.2748 


20.2843 


76 


2s 2 2p 2 3d 2 P 1/2 




22.0850 


22.0849 


22.0813 


35 




20.3736 


20.3735 


20.3840 


77 


2s 2 2p 2 3d 2 P 3/2 




22.1096 


22.1095 


22.1060 


36 




20.3755 


20.3754 


20.3831 


78 


2s 2 2p 2 3d 2 Si /2 




22.1807 


22.1806 


22.1599 


37 


2s2p33s 6 S 5/ 2 
2s 2 2p 2 3p 2 F 5/2 




20.6098 




20.7326 


79 


2s2p 3 3s 4 D 3/2 




22.3358 






38 




20.6772 


20.6771 


20.6810 


80 


2s2p 3 3s 4 D 5/2 




22.3365 






39 


2s 2 2p 2 3p 2 F 7/2 
2s 2 2p 2 3p 2 D 3/2 
2s 2 2p 2 3p 2 D 5/2 




20.6955 


20.6954 


20.7004 


81 


2s2p 3 3s 4 D 1/2 




22.3368 






40 




20.8957 


20.8956 


20.8739 


82 


2s2p 3 3s 4 D 7/2 




22.3429 






41 




20.8991 


20.8990 


20.8742 


83 


2s 2 2p 2 3d 2 D 5/2 




22.4934 


22.4933 


22.5332 


42 


2s 2 2p 2 3p 2 P 1/2 




21.0228 


21.0227 


21.0023 


84 


2s 2 2p 2 3d 2 D 3/1 




22.5124 


22.5123 


22.5509 



" Refers to the work of Bhatia & Landi |2003b). 



* Refers to data from the MCHF collection^!. 



For S 1 1+ , the present AS results agree well (within 20%) with 
predictions from other sources including superstructure (the 

data in chianti database, 94% of available transitions), MCHF^l 
(87%), MBPT (Merkelis et al. 199^1 83%) and the relativis- 
tic coupled-cluster theory (Nataraj et al. l2007h . for transitions 
between levels of the n-2 complex. For transitions from higher 
excited levels, e.g. n-3 configurations, only the unpublished cal- 
culation of Sampson & Zhang is available (from the chianti 
database.) Fig.[TJillustrates that only 46% of available transitions 
show agreement to within 20%. 

Additionally, we explicitly label some transitions with large 
differences in Fig.Q] They are all from the 50 th and 52 nd levels. 
We recall that we take configuration, total angular momentum 
and energy ordering to be the 'good' quantum numbers when 



4 http://open.adas.ac.uk 



level matching for comparisons. Exchanging the level matching 
for these two levels cannot eliminate the large difference now, 
unlike the case of S 9+ . Level mixing (2s2p3s 2 P contributes 86% 
for the 50 th , 2s2p3d 4 P contributes 90% for the 52 nd ) also can 
not explain this discrepancy for these strong transitions. 

Thus, we believe the atomic structure of the four iso-nuclear 
ions to be reliable, and expect the uncertainty in collision 
strengths due to in accuracies in the target structure to be cor- 
respondingly small. 

4. Scattering 

The scattering calculations were performed using a suite of 
parallel intermediate-coupling frame transformation (ICFT) R- 
matrix codes (Griffin et al. 1 1998b . We employed 40 contin- 
uum basis orbitals per angular momentum so as to represent 
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Table 5. Level energies (Ryd) of S ll)+ from different calculations along with experimentally derived values from NIST v4. Note: 
only levels with available NIST data and other predictions are listed. The complete table can be available from CDS archives. 



ID 


Specification 


3T1 

NISlLU 


AS 


LB03" 


MCHF'' 


ID 


Specification 


NIST 


AS 


LB03 


1 


2s 2 2p 2 J P 











48 


2s2p 2 3s 


J Po 




23.5543 




2 


2s 2 2p 2 3 Pj 


0.047459 


0.0468 


0.0493 


0.0463 


49 


2s2p 2 3s 


3 Pi 




23.5855 




3 


2s 2 2p 2 3 P 2 


0.112889 


0.1161 


0.1171 


0.1119 


50 


2s2p 2 3s 


3 P, 


23.6090 


23.6442 




4 


2s 2 2p 2 ! D 2 


0.611882 


0.6406 


0.6467 


0.6133 


51 


2s 2 2p3d 


'Pi 


23.5974 


23.7660 


23.6057 


5 


2s 2 2p 2 'S 


1.21134 


1.2242 


1.2376 


1.2114 


52 


2s 2 2p3d 


'F 3 


23.5958 


23.7845 


23.6205 


6 


2s2p 3 5 S 2 


1.69724 


1.6307 


1.6534 


1.6977 


56 


2s2p 2 3p 


5 D 2 


23.7613 


23.8851 




7 


2s2p 3 3 D 2 


3.23569 


3.2607 


3.2886 


3.2412 


64 


2s2p 2 3p 


3 D 3 


24.2379 


24.3250 




8 


2s2p 3 3 D! 


3.23832 


3.2616 


3.2911 


3.2439 


71 


2s2p 2 3s 


3 D 3 


24.5154 


24.5069 




9 


2s2p 3 3 D 3 


3.23819 


3.2686 


3.2920 


3.2438 


83 


2s2p 2 3d 


5 P 3 


25.0339 


25.1200 




10 


2s2p 3 3 P 


3.79950 


3.8307 


3.8625 


3.8049 


84 


2s2p 2 3d 


S P 2 


25.084 


25.1391 




11 


2s2p 3 3 Pi 


3.79986 


3.8335 


3.8632 


3.8052 


86 


2s2p 2 3d 


5 Pi 

L 1 


25.084 


25.1544 




12 


2s2p 3 3 P 2 


3.80380 


3.8391 


3.8676 


3.8089 


89 


2s2p 2 3d 


3 F 2 


25.174 


25.3020 




13 


2s2p 3 'D 2 


4.83133 


4.9362 


4.9602 


4.8413 


91 


2s2p 2 3d 


3 Fi 


25.207 


25.3358 




14 


2s2p 3 3 Si 


4.87728 


4.9635 


4.9810 


4.8859 


92 


2s2p 2 3d 


3 Fa 


25.257 


25.3825 




15 


2s2p 3 'Pi 


5.39908 


5.5099 


5.5368 


5.4087 


95 


2s2p 2 3p 


'D, 


25.420 


25.4352 




16 


2p 4 3 P 2 


7.39677 


7.4931 


7.5395 


7.4098 


97 


2s2p 2 3p 


'F 3 


25.470 


25.5594 




17 


2p 4 3 P, 


7.47723 


7.5708 


7.6232 


7.4897 


106 


2s2p 2 3d 


3 D 3 


25.514 


25.6540 




18 


2p 4 3 Po 


7.50561 


7.6000 


7.6530 


7.5193 


123 


2s2p 2 3d 


3 F 4 


26.266 


26.4556 




19 


2p 4 'D 2 


7.91401 


8.0543 


8.1063 


7.9285 


129 


2s2p 2 3d 


l F, 


26.353 


26.5317 




20 


2p 4 'So 


8.99180 


9.1670 


9.2198 


9.0064 


140 


2s2p 2 3d 


'D, 


26.6167 


26.8086 




21 


2s 2 2p3s 3 P 


21.2144 


21.2187 


21.0931 


21.1163 


158 


2s2p 2 3d 


3 P 2 


27.8214 


27.7608 




22 


2s 2 2p3s 3 Pi 


21.1438 


21.2440 


21.1182 


21.1405 


187 


2s 2 2p4s 


3 Pn 




28.7144 


28.5491 


23 


2s 2 2p3s 3 P 2 


21.2447 


21.3275 


21.2098 


21.2337 


189 


2s 2 2p4s 


3 Pl 




28.7282 


28.5618 


24 


2s 2 2p3s 'Pi 


21.3698 


21.4953 


21.3574 


21.3669 


191 


2s 2 2p4s 


3 P 2 




28.8229 


28.6645 


25 


2s 2 2p3p 'Pi 
2s 2 2p3p 3 D, 
2s 2 2p3p 3 D 2 




22.0422 


21.9292 


21.9462 


193 


2s 2 2p4s 


'Pi 




28.8753 


28.7036 


26 




22.1327 


22.0180 


22.0376 


197 


2s 2 2p4p 


3 Di 




29.0558 


28.8900 


27 




22.1520 


22.0356 


22.0560 


198 


2s 2 2p4p 


'Pi 




29.1039 


28.9380 


28 


2s 2 2p3p 3 D 3 




22.2253 


22.1128 


22.1348 


199 


2s 2 2p4p 


3 D 2 




29.1062 


28.9409 


29 


2s 2 2p3p 3 S, 
2s 2 2p3p 3 P 




22.2948 


22.1739 


22.1971 


200 


2s 2 2p4p 


3 Po 




29.1524 


28.9945 


30 




22.3509 


22.2032 


22.2262 


202 


2s 2 2p4p 


3 D 3 




29.1828 


29.0246 


31 


2s 2 2p3p 3 P, 




22.3907 


22.2521 


22.2756 


204 


2s 2 2p4p 


3 Si 




29.1886 


29.0242 


32 


2s 2 2p3p 3 P 2 
2s 2 2p3p 'D 2 




22.4250 


22.2820 


22.3050 


206 


2s 2 2p4p 


3 Pi 




29.2371 


29.0787 


33 




22.7280 


22.5690 


22.5600 


207 


2s 2 2p4p 


3 P 2 




29.2490 


29.0888 


34 


2s2p 2 3s 5 Pi 




22.9683 






208 


2s 2 2p4p 


'D 2 




29.3129 


29.1550 


35 


2s2p 2 3s 5 P 2 


22.9563 


23.0073 






211 


2s 2 2p4d 


3 F 2 




29.4298 


29.2586 


36 


2s2p 2 3s 5 P 3 


23.0130 


23.0624 






213 


2s 2 2p4p 


'So 




29.4357 


29.2712 


37 


2s 2 2p3p 'S 




23.0661 


22.8915 


22.8712 


216 


2s 2 2p4d 


3 F 3 




29.4670 


29.2987 


38 


2s 2 2p3d 3 F 2 




23.1255 


22.9832 


22.9897 


218 


2s 2 2p4d 


3 D 2 




29.4849 


29.3982 


39 


2s 2 2p3d 3 F 3 




23.1782 


23.0382 


23.0537 


221 


2s 2 2p4d 


3 D, 




29.5203 


29.3499 


40 


2s 2 2p3d 'D 2 


23.0757 


23.2087 


23.0685 


23.0775 


223 


2s 2 2p4d 


3 F 4 




29.5437 


29.3744 


41 


2s 2 2p3d 3 F 4 




23.2379 


23.1023 


23.1286 


226 


2s 2 2p4d 


'D 2 




29.5620 


29.3118 


42 


2s 2 2p3d 3 D, 


23.2229 


23.3724 


23.2202 


23.2203 


229 


2s 2 2p4d 


3 D 3 


29.458 


29.5925 


29.4288 


43 


2s 2 2p3d 3 D 2 


23.2349 


23.3858 


23.2360 


23.2428 


231 


2s 2 2p4d 


3 P 2 




29.6112 


29.4478 


44 


2s 2 2p3d 3 D 3 


23.2868 


23.4293 


23.2811 


23.2901 


232 


2s 2 2p4d 


3 P, 




29.6166 


29.4536 


45 


2s 2 2p3d 3 P 2 


23.3358 


23.4680 


23.3230 


23.3326 


233 


2s 2 2p4d 


3 Po 




29.6196 


29.4571 


46 


2s 2 2p3d 3 Pi 


23.3476 


23.4794 


23.3366 


23.3424 


238 


2s 2 2p4d 


'P, 




29.7068 


29.5459 


47 


2s 2 2p3d 3 P 




23.4854 


23.3448 


23.3505 


241 


2s 2 2p4d 


% 


29.5565 


29.7133 


29.5492 



MCHF 



23.6054 
23.6002 



" Refers to the work of Landi & Bhatia d2003h . 



the (N + 1) scattering electron for the four ions. All partial 
waves from 7=0 to 7=41 (S 9+ and S 11+ ) or 7=1/2 to 7=81/2 
(S 8+ and S 10+ ) were included explicitly and the contribution 
from higher 7-values were included using a "top-up" procedure 
(Burgess fl974l Badnell & Griffin l200ll) . The contributions from 
partial waves up to 7=12 (S 9+ and S 11+ ) or 7=23/2 (S 8+ and 
S 10+ ) were included in the exchange /^-matrix, while those from 
7=13 to 41 or 7=25/2 to 81/2 were included via a non-exchange 
/^-matrix calculation. In the exchange calculation, a fine energy 
mesh (1.0xl(T 5 z 2 Ryd, where z is the residual charge of ions) 
was used to resolve the majority of narrow resonances below the 
highest excitation threshold. From just above the highest thresh- 



old to a maximum energy of eight times the ionization potential 
for each ion, a coarse energy mesh (1.0xl0~ 3 z 2 Ryd) was em- 
ployed. For the non-exchange calculation, a step of 1.0xl0~ 3 z 2 
Ryd was used over the entire energy range. Additionally, ex- 
perimentally determined energies or adjusted energies were em- 
ployed in the MQDT expressions used by the ICFT method to 
further improve t he acc uracy of the results, as was done for 
Si 9+ (Liang et al., 2009a). The correction procedure was mainly 
done for levels of the n - 2 complex (needed because of the 
difficulty in obtaining a good structure here at the same time as 
describing n = 3 and 4 configurations with a unique orbital ba- 
sis) and some levels of the 2s 2 2p A 3s (where x=3,2,l or for 
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gf (present) gf (present) 



Fig. 1. Comparison of weighted oscillator stren gths gf of ele ctric-dipole transitions (to the 5 lowest-lying levels) for S 8+ - S 11+ . 
BL03 r efers to the work of Bh atia & Landi {2003a; 2003b), whereas LB03 and NSD07 correspond to the work of Landi & 
Bhatia (2003) and Nataraj et al. (120071) . respectively. The dashed lines correspond to agreement within 20%. [Colour online] 



S 8+ ' 9+ ' 10+ ' 11+ , respectively) configuration, as explained in detail 
in the structure section. 

We make use of the infinite energy Born limits (non-dipole 
allowed) or line strengths (dipole) to extend the /^-matrix col- 
lision strengths to higher scattering energies by interpo latio n 
of reduced variables, as described by Burgess & Tully (1992). 
Finally, thermally averaged collision strengths (T) were gener- 
ated at 13 electron temperatures ranging from 2xl0 2 (z + l) 2 K to 
2x 10 6 (z+ l ) 2 K. The data were stored in the ADAS adf04 format 
(Summers. |2004|) being available electronically from the OPEN- 
ADAS databaseEl, APAP-network0 and the CDS archives 

5. Results and Discussions 

5.1. S IX 

In Fig. [2] we make an extensive comparison of the present 
effecti ve colli sion strengths with the DW data of Bhatia & 
Landi (l2003al) for excitations from the ground level 2s 2 2p 4 3 P2. 
At the low temperature (log!T e (K)=5.1), only 27% of transi- 
tions show agreement within 20%. This can be easily explained 
by the o mission of resonances in the DW calculation by Bhatia 
& Landi d2003al) . At the temperature (logr e (K)=6.1) of peak 

5 http://www.apap-network.org 

6 http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 



fractional abundance in ionization equilibrium, the percentage 
is still low (41%). At the high temperature logr e (K)=7.1), the 
percentage increases to 58%. This is due to the reduced contribu- 
tions of near threshold resonances with increasing temperature. 
However, we note that there are a few transitions showing a ratio 
Tblib/Ticft >1.3, and the ratio increases with increasing tem- 
perature, e.g. the dipole transition of 2s 2 2p 4 3 P2 <— 2s 2 2p 3 3d 3 P; 
(1-80) marked by the dotted box in the lower-panel of Fig. [2] In 
the upper-panel of Fig. [2] we show the scaled collision strength 
as a function of reduced energy so as to shed light o n this o dd 
behaviour. The DW calculation by Bhatia & Landi (2003a|) is 
higher than the background of the present ICFT /^-matrix calcu- 
lation and the present Breit-Pauli DW (hereafter AS-DW) cal- 
culation using autostructure (Badnell 201 1). And the three dif- 
ferent calculations show a self-consistent behaviour approach- 
ing the infinite-energy limit point. So the odd behaviour is due 
to the higher background in the DW calculation by Bhatia & 
Landi d2003al) . The limit value from chianti (v6) is also plotted, 
which shows an excellent agreement with present calculations. 
This inconsistency in the chianti (v6) database is due to different 
data sources being adopted, e.g. the structure data is from a 24 
configuration calculation, whereas the scattering data is from a 
6 configuration calculation^. 



7 Landi, private communication (201 1). 
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Fig. 2. Comparison of (effective) collision strengths (Ts, see the lower panel) Q.s from the ground state of S 8+ . BL03 refers 
to the distorted wave calculation by Bhatia & Landi (2003a), AS-DW refers to the present Breit-Pauli DW calculation using 
autostructure. Upper panel: scaled collision strength for a dipole transition 2s 2 2p 4 3P2 <— 2s 2 2p 3 3d 3 Pi (1-80) with C=2.0. The 
limit value is 4g, /j / /£', ; - at 1.0 for the dipole transition. Lower panel: the ratio of Ts between the results of theDW calculation by 
Bhatia & Landi (2003a) and the ICFT /^-matrix calculation at log!T e (K) =5.1, 6.1 (corresponding to peak abundance of S 8+ in 
ionization equilibrium) and 7.1. The dashed lines correspond to agreement within 20%. The transition marked by dotted box is the 
dipole transition 1-80 shown in the upper panel. [Colour online] 
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Fig. 3. Comparison of effective collision stre ngths of S 8+ with 
the /^-matrix results of Butler and Zeippen (1 19941) for transi- 
tions of the ground configuration 2s 2 2p 4 . Fi lled sy mbols with 
solid curves are results of Butler and Zeippen (1 19941) . while open 
symbols with dotted curves corresponds to the present ICFT R- 
matrix calculation. Note: The same symbol in the two sets of 
results corresponds to the same transition. [Colour online] 



For S 8+ , an earlier /^-matrix calculation for transitions within 
the ground configuration is available (Butler & Zeippen, 1994) 
for which the LS-coupling K-matrices were transformed alge- 



braically to intermediate coupling to obtain collision strengths 
between the fine- structure levels. A detailed comparison has 
been made between the two different /^-matrix calculations, see 
Fig. [3] At the low temperature (r e ~ 1.0 x 10 4 K), there is 
a large difference between the two different /^-matrix calcula- 
tions. A separate ICFT 7?-matrix calculation with finer mesh 
(1.0x1 0~ 6 z 2 ) near threshold confirms that the effect of resonance 
resolution is less than 2% for nearly all excitations, except for 
the 2-5 (10% at logr e (K)=4.1) and 3-5 (24% at logr e (K)=4.1) 
transitions. So the present effective collision strengths are gen- 
erally converged with respect to resonance resolution. The large 
differences between the two different /^-matrix calculations may 
be due to deficiencies in the transformational approach us ed by 
Butler & Zeippen d 19941) . as deta iled by Griffin et al. (1998) 
and demonstrated by Liang et al. d2008l) . The adoption of ob- 
served energies for levels of n-2 complex in the present ICFT R- 
matrix calculation gives better positioning of near threshold res- 
onances tha n the p revious ones with theoretical energies (Butler 
& Zeippen, |1994|) . So, the present effective collision strengths 
are expected to be more reliable at low temperatures. 



5.2. S X 

In Fig. [4] an extensive comparison has been made with previ- 
ous ^-matrix calculation (Bell & Ramsbottom 2000) at three 
temperatures: logr e (K)=5.2, 6.2 (corresponding to peak frac- 
tion in ionization equilibrium) and 6.7. At the low temper- 
ature (logr e (K)=5.2), only 27% of all available transitions 
show an agreement within 20%. Even at the high tempera- 



8 



Liang et al.: ^-matrix calculation of four iso-nuclear sulphur ions 



Table 7. Comparison of weighted oscillator streng ths (g f) of 
S 9+ between the previous data (Bhatia & Landi l2003al) and 



the present autostructure calculations with/without valence- 
valence two-body fine-structure interactions (TBFS) for the 
ground complex and level energy correction (labeled as LEC). 
The index number corresponds to that in Table [4] The last col- 
umn is the data presented in Fig [1] Note: only data with differ- 
ence being > 20% are listed. 

Present 



With TBFS 



BL03 



LEC 



Without TBFS 
LEC 



10 
13 
19 
20 
32 
56 
76 
77 
83 
84 
7 
20 
59 
7 
16 
17 
8 
9 
17 
18 
51 
58 
59 
61 



4.176-6 
1.195-4 
1.214-4 
1.526-4 
2.007-6 
1.386-8 
8.706-5 
7.970-4 
1.856-7 
3.191-7 
1.198-5 
8.382-4 
2.976-5 
2.021-6 
4.998-5 
1.237-4 
1.704-6 
1.227-4 
5.655-4 
2.940-6 
2.586-5 
1.141-3 
3.933-4 
3.205-4 



4.240-6 
1.222-4 
1.207-4 
1.518-4 
1.979-6 
1.478-8 
8.684-5 
7.942-4 
1.603-7 
3.229-7 
1.196-5 
8.311-4 
2.927-5 
2.014-6 
4.956-5 
1.222-4 
1.686-6 
1.256-4 
5.599-4 
2.848-6 
2.574-5 
1.133-3 
3.888-4 
3.171-4 



1.330-5 
1.410-4 
1.483-4 
1.951-4 
1.035-6 
9.315-7 
6.689-5 
6.280-4 
1.601-5 
5.315-7 
1.608-5 
6.546-4 
3.587-5 
2.426-6 



6.865- 

1.487- 

9.851- 

1.323- 

6.878- 

2.795- 

1.494-5 

8.269-4 

5.672-4 

4.305-4 



6.685-6 
1.392-4 
1.232-4 
1.415-4 
1.315-6 
3.075-7 
8.707-5 
7.364-4 
2.091-6 
3.216-0 
1.347-5 
8.107-4 
2.142-5 
2.116-6 
4.955-5 
1.222-4 
1.309-6 
1.264-4 
5.443-4 
4.540-6 
2.086-5 
1.026-3 
4.425-4 
3.524-4 



9.749-6 
1.225-4 
1.455-4 
2.084-4 
1.722-6 
2.317-7 
6.677-5 
6.882-4 
7.996-6 
1.255-8 
1.464-5 
7.251-4 
4.579-5 
2.311-6 
6.868-5 
1.486-4 
1.323-6 
1.625-4 
7.040-4 
1.326-6 
2.190-5 
9.380-4 
5.000-4 
3.855-4 



Note: x±y = xx I0 ±y . 



ture (logr e (K)=6.7), the percentage is only about 34%. Ratios 
(t^t) less than unity can be understood in terms of the 
finer energy mesh used (present: 1.0xl0~ 5 z 2 Ryd, Bell & 
Ramsbottom l2000t >0.008 Ryd) and resonances attached to the 
2s 2 2p 2 3/ configurations in our present ICFT R-matiix calcula- 
tion, as well as the purely algebraic jajom approach that was used 
by Bell & Ramsbottom (2000). However, the ratio being larger 
than unity requires another explanation. So, we select one tran- 
sition marked by the bold '\' in Fig. [4] to investigate the source 
of the difference between the two different /^-matrix results. 

Figure [5] shows a comparison of our present (effective) 
collision strength with the previous 7?-matrix results for the 
2s 2 2p 3 4 S3/2 <— 2s2p 4 2 P3/2 (1-12) dipole transition. Around 
the temperat ure of T e ~ 9.0 x 10 5 K - 6.0xl0 6 K, the Bell & 
Ramsbottom (2000) result is higher than present ICFT /^-matrix 
calculation, and by up to a factor of 2. We note that some 
pseudo-orbitals (3p, 4s, 4d and 4f) were included in the work 
of Bell & Ramsbottom (2000). They stated that some pseudo- 
resonances are found above the highest threshold (19.682 Ryd). 
One of the authors (Ramsbottom, private communication, 201 1) 
has provided us with collision strengths (Q) with the pseudo- 
resonances at high energies removed. A comparison of the scaled 
collision strengths Q reveals that the backgrounds of the two 
different /^-matrix calculations agree well, an d are consistent 
with the DW calculation by Bhatia & Landi (l2003bl) . So, the 
large difference between the two different /^-matrix calculations 
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Fig. 4. Comparison of effective collision strengths of S 9+ with 
the jajom /^-matrix results of Bell & Ramsbottom (2000J for 
transitions of the n = 2 complex and 2s 2 2p 2 3s configuration. 
One point marked by bold '\' refers to the 2s 2 2p 3 4 S3/2 <— 
2s2p 4 2 P3/2 transition (1-12), that will be examined in Fig. [5] 
[Colour online] 
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Fig. 5. Comparison of excitation data for the 2s 2 2p 3 4 S3/2 
<- 2s2p 4 2 P 3/2 transition (1-12) of S 9+ . Here, BR00 refer s 
to the /^-matrix calculation of Bell & Ra msbotto m yOOO), 
BL03 to the DW result of Bhatia & Landi d2003bl) and ICFT 
to the present calculation. Upper panel: Effective collision 
strengths — the BR00 without pseudo-resonances result was re- 
derived by us from the said original collision strengths provided 
by Ramsbottom (private communication, 2011). Lower panel: 
Scaled collision strengths, with the scaling parameter C set to 
2.0. [Colour online] 



is not arising from the difference in their structures. We then 
re-derived the effective collision strengths, which shows the 
expected behaviour, see Fig. [5] So, it appears that the previ- 
ously published 7?-matrix effective collision strengths of Bell & 
Ramsbottom (2000) were derived from their collision strengths 
before the pseudo-resonances were subtracted. So, the ratios 
greater than unity in Fig.|4]should be mostly/partly attributed to 
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the pseudo-resonances in the previous /^-matrix calculation. So, 
the present results are more reliable for modelling applications. 

For excitations to higher levels of n=3 configurations, 
only DW data i s available, e.g. the latest work of Bhatia 
& Landi (20033). A comparison there demonstrates that the 
resonance contribution is strong for some transitions and is 
widespread, as expected. For conciseness, the figure is not shown 
here. 

5.3. S XI 

As mentioned in the introduction, interpolated data f rom R- 
matrix results is available for S 10+ (Conlon et al. Il992l) . These 
resultant data are valid over a temperature range approximately 
equal to T e ~ 3.2 x 10 5 — 1.3 x 10 7 K for S 10+ . In this temper- 
ature range, the interpolated excitation data show a good agree- 
ment with present the ICFT /^-matrix calculation for almost all 
transitions, as shown in Fig. [6] even though only partial waves 
of L <9 and 12 continuum basis orbitals in each channel were 
included. That is, the effective collision strength is converged 
in this temperature range using a small range of partial waves 
etc. Lennon & Burke (1994, hereafter LB94) performed an R- 
matrix calculation which included all 12 terms of the ground 
complex and adjusted the diagonal elements of the LS-coupling 
Hamiltonian matrix to the (fine-structure averaged) observed en- 
ergies before diagonalization. They provided data for transitions 
between fine-structure levels in the ground configuration 2s 2 2p 2 
plus the 2s2p 3 5 S 2 . At low temperatures T e < 1.0 x 10 5 K, the 
present ICFT /^-matrix calculation is systematically higher than 
this previous small-scale /^-matrix result except for the 1-6 tran- 
sition, see Fig. [6] This situation can likely be attributed to the 
much larger close-coupling expansion (to n — 4) and associated 
resonances in the present calculation. We recall also that we used 
observed level energies in the present ICFT /^-matrix calculation 
via multi-channel quantum defect theory (MQDT). 

In case of the 1 - 6 tran sition, the original collision strength 
of Lennon & Burke (1 19941) is available from TIPTOPbaseB In 
figure [7] we compare the two sets of results. We see that there is 
a somewhat oddly high background around 0.5-2.0 Ryd in the 
results of Lennon & Burke (1994). This is the likely reason their 
effective collision strength is notably larger than the present one 
at lower temperatures. 

Comparison with the DW calculation of Landi & 
Bhatia d2003l) demonstrates that only 22% of all available tran- 
sitions show agreement within 20%. Fig. [8] demonstrates that 
the resonance contribution is strong for some transitions, and 
widespread as expected again. At high temperature, uncertain- 
ties of scattering data are dominated by the accuracy of struc- 
ture calculation because the resonance contribution becomes 
increasingly small. But only 43% of all available transitions 
show agreement within 20%, which is significantly lower than 
that in the assessment for weighted oscillator in section [3] We 
also notice there are a few transitions showing the ratio be- 
ing larger than unity. So we select one transition (2s 2 2p 2 3 Po 
<— 2s 2 2p4p 3 Po, see the bold '\' mark in Fig. [8}a) to inves- 
tigate. Fig. [8] clearly demonstrates that the DW data of Landi 
& Bhatia (2003) is higher than the background of the present 
ICFT /?-matrix calculation. But the pr esent B reit-Pauli DW cal- 
culation Using AUTO STRUCTURE (BadnellHon]) shows an excellent 
agreement with the background of the /^-matrix calculation — 
both us e the e xact same atomic structure. As stated by Landi & 
Bhatia d2003h . a small atomic model (nine lowest configurations, 
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Fig. 7. Collision strengths (Q) of the 2s 2 2p 2 3 P <- 2s2p 3 5 S 2 
(1-6) transitions of S 10+ , where LB94 corresponds to the R- 

matrix work of Lennon & Burke (11994 from the TIPTOPbaseP, 
ICFT denotes the present work. [Colour online] 

72 fine-structure levels) was adopted in their scattering calcula- 
tion because of their available computer resource. So, we per- 
formed another separate AS-DW calculation with the 9 lowest 
configurations, in which the optimization procedure is done as 
mentioned above for S I0+ . The resultant data show good agree- 
ment with the DW calculation by Landi & Bhatia d2003l) . So 
ratios lower than unity and the low percentage of agreement in 
the scatter plot mentioned above are likely due to the use of a 
much larger configuration interaction expansion in the present 
ICFT /^-matrix calculation. 

5.4. S XII 

As stated by Keenan et al. d2002l). a small error in the previous 
excitation data (Zhang et al. 1994) was found for a few transi- 
tions of some boron-like ions, and those data were replaced. In 
Fig |9j we compare the present ICFT ^-matrix excitation data 
with the revised data of Keenan et al. ( 2002) at three different 
temperatures (logr e (K)=6.04, 6.40 and 6.78) to check the va- 
lidity of the present results or improvement by including larger 
CI and extensive close-coupling expansions. For strong exci- 
tations (>0.1), a good agreement (within 20%) is obtained for 
most excitations (82%). For weak excitations, the present ICFT 
/^-matrix results are systematically larger than previous ones ex- 
cept for a few transitions, e.g. 8-13 and 9-13. Indeed, the weaker 
the excitation, the greater the difference, and by more than a 
factor of 2 for a group of the weakest excitations. This can be 
easily explained by resonances attached to «=3 levels included 
in the present work, and this effect is stronger for weaker ex- 
citations. For the two above mentioned transitions (8-13 and 
9-13), the previous /^-matrix calculation is significantly higher 
than the present ones at logr e (K)=6.04 by a factor of 2.5 and 
40%, respectively. Unfortunately, there are no previous colli- 
sion strengths available to compare with — examination of the 
present collision strengths uncovers no untoward behaviour for 
these two transitions. 

For excitations to higher excited levels of the n=3 configura- 
tions, only an unpublished DW calculation (Zhang & Sampson 
1995) is available — compiled in the chianti database. A com- 
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parison demonstrates that the resonance contribution is strong 
for some transitions, and is widespread as expected. For con- 
ciseness, the figure is not shown here. 



Additionally, we checked the sensitivity of the high-T e Ts 
to the top-up and find that it is greatest on the weakest (dipole) 
transitions but it is not significant compared to the inherent un- 
certainties in the atomic structure (f-values) for such transitions 
— the strong transitions are well converged. 



6. Summary 

Electron-impact excitation data for four iso-nuclear sulphur ions 
(S 8+ , S 9+ , S 10+ and S 1I+ ) have been calculated using the ICFT 
/^-matrix method with extensive CI and large close-coupling ex- 
pansions, as listed in Table [TJ 

Good agreement overall with the available experimentally 
derived data and other theoretical results for level energies 
and weighted oscillator strengths supports the reliability of the 
present /{-matrix excitation data. 

For excitations to levels of the n-2 complex, an extensive 
assessment have been made with previous /^-matrix calculations 
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Fig. 9. Comparison of effective collision strengths of S 11+ for 
all excitations between levels of the n-2 complex at three differ- 
ent temperatures log7' e (K)=6.04, 6.40 and 6.87. Ticft refers to 
the present ICFT /^-matrix calculation, and KKR0 2 corr esponds 
to the previous /^-matrix results by Keenan et al. (|2002|) . A few 
transitions with large difference are marked by labels around the 
points, and are linked together for results at different tempera- 
tures. [Colour online] 
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Keenan, F.P, O'Shea, E., Thomas, R.J., et al. 2000, MNRAS, 315, 450 
Landi, E., & Bhatia, A.K. 2003, ApJS, 149, 251 
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Liang, G.Y, Whiteford, A.D., & Badnell, N.R. 2008, J. Phys. B: At. Mol. Opt. 
Phys., 41, 235203 

Liang, G.Y, Whiteford, A.D., & Badnell, N.R. 2009a, A&A, 499, 943 
Liang, G.Y, Whiteford, A.D., & Badnell, N.R. 2009b, A&A, 500, 1263 
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51,233 

Nataraj, H.S., Sahoo, B.K., Das, B.P, et al. 2007, J. Phys. B: At. Mol. Opt. 
Phys., 40, 3153 

Raassen, A.J.J., Mewe, R., Audard, M., et al. 2002, A&A, 389, 228 
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Tachiev, G.I., & Froese Fischer, C. 2002, A&A, 385, 716 
Thomas, R.J., & Neupert, W.M. 1994, ApJS, 91, 461 

Witthoeft, M.C., Whiteford, A.D., & Badnell, N.R. 2007, J. Phys. B: At. Mol. 
Opt. Phys., 40, 2969 
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available to check the validity and improvement of the present 
ICFT /^-matrix results. For excitations to higher excited levels of 
n=3 and/or 4 configurations, only DW calculations are available 
to compare with. The improvement of the present calculations is 
illustrated as expected by including resonances. For some tran- 
sitions, configuration interaction has a significant effect on the 
atomic structure and this carries through to the final (effective) 
collision strengths, as shown in the cases of S 8+ and S 1()+ . 

In conclusion, the present ICFT /^-matrix excitation data of 
S>8+,9+,io+ andS 11+ are assessed to be valid over an extensive tem- 
perature range, and a significant improvement is achieved over 
previous available ones to date due to the extensive CI and large 
close-coupling expansions used in the present work. This will 
replace data from DW and small /^-matrix calculations presently 
used by astrophysical and fusion communities, and its use can be 
expected to identify new lines, improve spectral analyses and di- 
agnostics of hot emitters or absorbers in astrophysics and fusion 
researches. 
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Table 6. Level energies (Ryd) of S 11+ from different calculations along with experimentally derived values from NIST v4. Note: 
only levels with available NIST data and other predictions are listed. The complete table can be available from CDS archives. 



ID 


Specification 


zm — 


AS 


CHIANTP 


NSD07 6 


ID 


Specification 


NIST 


AS 


CHIANTI" NSD07* 


1 


2s 2 2p 2 P 1/2 
2s 2 2p 2 P 3/2 




0.0000 






54 


2s2p3d 2 F 5/2 


27.322 


27.4185 


27.2824 


2 


0.119698 


0.1185 


0.1 197* 


0.1227 


55 


2s2p3d 2 F 7/2 


27.396 


27.4896 


27.3578 


3 


2s2p 2 4 P I/2 
2s2p 2 4 P 3/2 


1.76678 


1.7295 


1.7668' 


1.7656 v 


56 


2s2p3d 2 P 3/2 


27.443 


27.5264 


27.3679 


4 


1.81046 


1.7733 


1.8105* 


1.8093 v 


57 


2s2p3d 2 P l/2 


27.477 


27.5669 


27.4098 


5 


2s2p 2 4 P 5/2 


1.87197 


1.8394 


1.8720* 


1.8702* 


58 


2s2p3p 2 P 1/2 




27.8631 


27.7804 


6 


2s2p 2 2 D 3/2 


3.1594 


3.1941 


3.1594* 


3.146P 


59 


2s2p3p 2 P 3/2 




27.8800 


27.7947 


7 


2s2p 2 2 D 5/2 


3.16214 


3.1988 


3.1621* 


3.1493 v 


60 


2s2p3p 2 D 3/2 


27.894 


27.9040 


27.8174 


8 


2s2p 2 2 S 1/2 
2s2p 2 2 P 1/2 
2s2p 2 2 P 3/2 


4.0057 


4.0522 


4.0057* 


3.9887 v 


61 


2s2p3p 2 D 5/2 


27.894 


27.9057 


27.8174 


9 


4.23516 


4.3012 


4.2352* 


4.2260* 


62 


2s2p3p 2 Si /2 




28.1199 


28.0381 


10 


4.2960 


4.3666 


4.2960* 


4.2864 v 


63 


2p23s 4 P 1/2 




28.3086 


28.2094 


11 


2p 3 4 S 3/2 


5.55941 


5.5632 


5.5594* 


5.5529>- 


64 


2 P 23s 4 P 3/2 




28.3510 


28.2527 


12 


2p 3 2 D 3/2 


6.2869 


6.3534 


6.2869* 


6.2603 v 


65 


2p23s 4 P 5/2 




28.4161 


28.3157 


13 


2p 3 2 D 5/2 


6.2921 


6.3634 


6.2921* 


6. 2645 y 


66 


2s2p3d 2 F 7/2 




28.6562 


28.5654 


14 


2p 3 2 Pi /2 


7.0534 


7.1496 


7.0534* 


7.0243 v 


67 


2s2p3d 2 F5/ 2 


28.506 


28.6584 


28.5682 


15 


2p 3 2 P 3/2 


7.06966 


7.1678 


7.0697* 


7.0398 v 


68 


2s2p3d 2 D 3/2 




28.7631 


28.6427 


16 


2s23s 2 Si /2 




23.2974 


23.1956 


23.3186 


69 


2s2p3d 2 D 5/2 


28.620 


28.7763 


28.6548 


17 


2s23p 2 Pi /2 




24.1824 


24.0825 


24.2053 


70 


2p23s 2 P 1/2 




28.8013 


28.6888 


18 


2s23p 2 P 3/2 




24.2153 


24.1162 


24.2386 


71 


2p23s 2 P 3/2 




28.8704 


28.7571 


19 


2s23d 2 D 3/2 


25.036 


25.0526 


24.9441 


25.0837 


72 


2s2p3d 2 P l/2 




28.9214 


28.8042 


20 


2s23d 2 D 5/2 


25.043 


25.0624 


24.9536 


25.0926 


73 


2s2p3d 2 P 3/2 




28.9346 


28.8142 


21 


2s2p3s 4 Pi/ 2 




25.1340 


25.0355 




74 


2p23p 2 S 1/2 




28.9545 


28.8417 


22 


2s2p3s 4 P 3 / 2 




25.1728 


25.0740 




75 


2p23p 4 D 1/2 




29.0458 


28.9315 


23 


2s2p3s 4 P5/ 2 


25.285 


25.2443 


25.1480 




76 


2p23p 4 D 3/2 




29.0716 


28.9572 


24 


2s2p3s 2 Pj/ 2 




25.5847 


25.4734 




77 


2p23s 2 D 5/2 




29.0868 


29.0118 


25 


2s2p3s 2 P 3 / 2 




25.6643 


25.5546 




78 


2p23s 2 D 3/2 




29.0915 


29.0156 


26 


2s2p3p 4 Di/ 2 




25.9258 


25.8279 




79 


2p23p 4 D 5/2 




29.1157 


29.0016 


27 


2s2p3p 4 D 3/2 




25.9617 


25.8625 




80 


2p23p 4 D 7/2 




29.1757 


29.0608 


28 


2s2p3p 2 Pi /2 


25.886 


25.9950 


25.9018 




81 


2 P 23p 4 P 1/2 




29.1961 


29.0807 


29 


2s2p3p 2 P 3/2 


25.953 


26.0102 


25.9171 




82 


2p23p 4 P 3/2 




29.2084 


29.0907 


30 


2s2p3p 4 D 5/2 




26.0225 


25.9227 




83 


2p23p 4 P 5/2 




29.2434 


29.1265 


31 


2s2p3p 4 D 7/2 




26.0858 


25.9892 




84 


2p23p 2 D 3/2 




29.3169 


29.1974 


32 


2s2p3p 4 S 3/2 




26.1845 


26.0807 




85 


2p23p 2 D 5/2 




29.3979 


29.2818 


33 


2s2p3p 4 Pi/ 2 




26.2653 


26.1647 




86 


2p23p 2 P 3/2 




29.5566 


29.4369 


34 


2s2p3p 4 P 3/2 




26.3044 


26.2043 




87 


2p23p 2 P 1/2 




29.5699 


29.4480 


35 


2s2p3p 4 P 5/2 




26.3375 


26.2349 




88 


2p23p 4 S 3/2 




29.6045 


29.3999 


36 


2s2p3p 2 D 3/2 


26.388 


26.4104 


26.3067 




89 


2p23d 4 F 3/2 




29.7270 


29.6128 


37 


2s2p3p 2 D 5/2 


26.468 


26.4831 


26.3827 




100 


2p23d 2 F 5/2 




30.0001 


29.8983 


38 


2s2p3d 4 F 3 / 2 




26.6901 


26.5796 




101 


2p23d 2 P 1/2 




30.0173 


29.8852 


39 


2s2p3d 4 F 5/2 




26.7135 


26.6030 




102 


2p23d 2 F 7/2 




30.0834 


29.9856 


40 


2s2p3p 2 Si/ 2 


26.698 


26.7294 


26.6255 




103 


2p23p 2 D 3/2 




30.1533 


29.9758 


41 


2s2p3d 4 F7/ 2 




26.7479 


26.6388 




104 


2p23p 2 D 5/2 




30.1562 


29.9741 


42 


2s2p3d 4 Fc>/ 2 




26.7970 


26.6920 




105 


2p23d 4 P 5/2 




30.1599 


30.0059 


43 


2s2p3d 4 D 1/2 




26.9396 


26.8016 




122 


2p23p 2 P 1/2 




31.1471 


31.0749 


44 


2s2p3d 4 D 3/2 


26.903 


26.9407 


26.8021 




123 


2p23p 2 P 3/2 




31.1607 


31.0883 


45 


2s2p3d 4 D 5/2 


26.890 


26.9417 


26.8024 




129 


2s24d 2 D 3/2 


32.289 


32.3629 




46 


2s2p3d 2 D 3/2 




26.9844 


26.8391 




130 


2s24d 2 D 5/2 


32.289 


32.3668 




47 


2s2p3d 4 D 7/2 


26.952 


26.9908 


26.8553 




149 


2s2p4p 2 D 5/2 


34.008 


34.1315 




48 


2s2p3s 2 P 1/2 




26.9965 


26.9334 




159 


2s2p4d 4 D 7/2 


34.209 


34.2968 




49 


2s2p3d 2 D 5/2 


26.943 


26.9981 


26.8553 




172 


2s2p4d 2 F 5/2 


34.325 


34.4282 




50 


2s2p3s 2 P 3 / 2 




27.0021 


26.9344 




177 


2s2p4d 2 F 7/2 


34.370 


34.4854 




51 


2s2p3d 4 P 5/2 


27.008 


27.0528 


26.9209 




193 


2s2p4d 2 F 7/2 


35.84 


35.9979 




52 


2s2p3d 4 P 3/2 


27.036 


27.0622 


26.9367 




194 


2s2p4d 2 F 5/2 


35.84 


35.9981 




53 


2s2p3d 4 P 1/2 




27.0702 


26.9457 




204 


2s2p4f 2 D 5/2 




36.2258 





8 Refers to the theoretical value in chianti v6. 

* NSD07 corresponds to the prediction of Nataraj et al. d2007h with relativistic coupled-cluster theory. 

* These data are found to be observed, after checking the original paper (Zhang et al. ll994h . 

1 Data are from the work of Merkelis e al. d!995t) by using many-body perturbation theory (MBPT). 
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